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ABSTRACT 

The early optical emission of the moderately high redshift (z = 3.08) GRB 060607A 
shows a remarkable broad and strong peak with a rapid rise and a relatively slow 
power-law decay. It is not coincident with the strong early-time flares seen in the X- 
ray and gamma-ray energy bands. There is weak evidence for variability superposed 
on this dominant component in several optical bands that can be related to flares in 
high energy bands. While for a small number of GRBs, well-sampled optical flares 
have been observed simultaneously with X-ray and gamma ray pulses, GRB 060607A 
is one of the few cases where the early optical emission shows no significant evidence 
for correlation with the prompt emission. In this work we first report in detail the 
broad band observations of this burst by Swift. Then by applying a simple model 
for the dynamics and the synchrotron radiation of a relativistic shock, we show that 
the dominant component of the early emissions in optical wavelengths has the same 
origin as the tail emission produced after the main gamma ray activity. The most 
plausible explanation for the peak in the optical light curve seems to be the cooling of 
the prompt after the main collisions, shifting the characteristic synchrotron frequency 
to the optical bands. The fact that the early emission in X-ray does not show a 
steep decay, like what is observed in many other GRBs, is further evidence for slow 
cooling of the prompt shell within this GRB. It seems that the cooling process requires 
a steepening of the electron energy distribution and/or a break in this distribution 
at high energies. From simultaneous gamma-ray emission during the first flare, the 
behaviour of hardness ratio, and the lack of spectral features, we conclude that the 
X-ray flares are due to the collision of late shells rather than late reprocessing of the 
central engine activities. The sharp break in the X-ray light curve at few thousands of 
seconds after the trigger, is not observed in the IR/optical/UV bands, and therefore 
can not be a jet break. Either the X-ray break is due to a change in the spectrum of 
the accelerated electrons or the lack of an optical break is due to the presence of a 
related delayed response component. 
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1 INTRODUCTION 

GRB 060607A was a long, fairly hard GR B localized by the 
Swift Burst Alert Telescope (BAT K|Barthelmv et ail 
l2005al fbh at T ~ 05:12:13 UT on 2006 June 
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7 l|Ziaeepour et all [2006). The Swift spac ecraft rapidly 
slewe d, directing the X-Ray Telescope (XRT) l|Burrows et al.l 
2005) and the Ult r aviolet and Optical Telescope 
fUVOT) (|Roming et al.l 120051 ) at the BAT position. 
The observations commenced at To + 63.6 sec after the BAT 
trigger in the optical/UV and To + 73.6 sec in the X-ray, 
providing broad-band spectral and high time resolution 
light curves in the X-ray, optical, and ultraviolet. The 
X-ray light curve shows two bright flares, the first of which 
at To + 98 sec is also observed in the BAT energy range 
15 — 300 keV. The second flare occurs at To + 260 sec and 
is marginally (if at all) detected by BAT (see Fig. HJ). The 
optical emission is at first quite faint, then rises nearly three 
magnitudes over a time span of about 100 sec. The peak in 
the optical is not simultaneous with the X-ray flares. 

Historically, ground-based follow-up optical observa- 
tions of GRBs usually started at least several minutes after 
the burst occurred, and it has not been unusual for optical 
observations to start several hours into the afterglow phase. 
Over the past few years fast-slewing robotic telescopes e.g . 
llBoerll200ll:|Perez- Ramirez et al.ll2"oo3;IVestrand et al.ll2002l : 
lAkerlof et all 120031; ICovino et all 1 20041 ) and the Swift mis- 
sion l|Gehrels et al.l 12004 ) have made it possible to observe 
the first few minutes o f optical emiss i on fro m GRBs. In some 
cases: GRB 990123 (lAkerlof et all Il999l ). GRB 041219A 
JVestrand et al l I2005I ), GRB 05 0820A (IVestrand et all 
l2006l ), GRb"051109A, GRB 051111 | Yost et al.ll2007l), GRB 
06012 4 (|Romano et al.1120061). GRB 0604 18 l|Molinari et all 
120071 ). GRB 060526 JPai et al l 120071 ), GRB 060607 A 
GRB 061007 i|Schadv et al.ll2006l), GRB 0611 21 (|Page et all 
l2007h . GRB 070616 (|Starling et al l l2007al lbh the optical 
emission has been observed during the prompt gamma-ray 
emission. In the case of GRB 990123 there is no correlation 
between the structure of the ga mma-ray light curve and the 
optical flux (jAkerlof et al.ll 19991 ) . If anything, there is an ap- 
parent anti-correlation, although the optical data is not well 
sampled. On the other hand, all the other bursts above show 
a correlation between the optical and gamma-ray fluxes. In 
the case of GRB 060124, although the optical flux slightly 
increases during the main peak at about 600 sec, its rela- 
tive rise is much smaller than the rise in gamma-ray and 
X-ray fluxes. In the case of GRB 060607A, clearly there is 
no correlation between the dominant component of the early 
optical and X-ray emissions. The optical flux arrived at its 
peak when the X-ray emission was decreasing. A similar be- 



havio ur was also observed in GRB 060418 (|Molinari et alj 
l2007i ). 

It is not yet clear how the prompt gamma-ray emis- 
sion is related to the later emissions in the lower en- 
ergy bands. Swift XRT observations have shown that the 
X-ray early emission of the most GRBs exhibit a very 
steep power-law decline with (3 < a < 5, where the flux 

at time t is f„ oc t~ a ) decay within Tq+ < 1000 sec 

llNousek et alj|2006l; ICTBrien et al.|[2006l ; iButler et al1l2006l; 
IWillingale et al.l 120071 )). This decay was expected to be 
caused by the high-latitude emission from the internal 
shocks that are driving the prompt high-energy emission, 
and therefore is not produced by t he same mechanism 
that drives the late-tim e afterglow (e.g. Zh ang fc Kobavashil 
l2005l ; lLiang et al.|[2006l ). presumably the external shock with 
the ISM or wind. However, this rapid decay phase has not 
been observed in any optical early emission. Moreover, the 



high-latitude emission should satisfy a str ict relation be- 
tween a and the photon index (3: a — 2 + /3 (Fen imore et al.l 
1996). This relation is however only satisfied in a small 
fraction of GRBs with a steeply decaying X-ray tail emis- 
sion which also do not show significant spectral evolution 
l|Zhang et al. 2007a). The spectral evolution of the tail emis- 
sion is another argument against high latitude origin of the 
tail. Therefore, one can conclude that the high-latitude emis- 
sion is not a dominant contributor in the tail emission. The 
fact that the light curves in different bands in general do not 
follow each other suggest s that either multiple components 
should be involved l|Kumar et al.|[2006l ) and/or processes are 
chromatic and the emission evolves both in time and in en- 
ergy. Therefore, one should expect various degrees of corre- 
lation between the light curves in the different energy bands 
depending on the internal properties of the system and its 
environment. GRB 060607A is a good example of how pe- 
culiar the relation between energy bands can be. The initial 
decay slope in X-ray is shallower than most bursts a ~ 1, 
and the optical light curve shows a peak uncorrelated to the 
flares in gamma-ray and X-ray. 

In this paper we focus on the early optical and X-ray 
emissions as observed by Swift's UVOT and XRT instru- 
ments. The goal is to compare the multi-band Swift data, 
with a simple relativistic shock and synchrotron emission 
model, and to try to reconstruct the history of events lead- 
ing to the prompt and afterglow emission as observed by 
Swift. In Section [2] we describe the Swift observations. In 
Section[3l we present the time variation analysis, and in Sec- 
tion rj] we discuss broad-band spectral variability during the 
early emission phase and give a qualitative interpretation of 
the data. In Section we apply the model mentioned above 
to the data to interpret observations and to estimate some of 
the parameters. Finally, we summarize our results in Section 
[6] The theoretical model used in Section[5]is reviewed briefly 
in Appendix IA1 In Appendix iBl the conditions determining 
the reliability of the extrapolation of BAT light curve to the 
XRT energy band is discussed. 



2 OBSERVATIONS 

In this section we briefly report the result of the analysis 
of the BAT, XRT, and UVOT data. Due to the peculiarity 
of the optical afterglow of this burst and its importance for 
the interpretation of the GRB, the UVOT data is discussed 
with more details. 



2.1 BAT Data 

The BAT light curve showed a double-peaked str ucture 
with a duration of about 40 sec (|Tueller et alj |2006). The 
peak count rate was approximately 3000 count sec -1 (15 — 
350 keV) at the time of the trigger. The mask- weighted light 
curve (Fig. [TJ consists of two overlapping FRED-like peaks 
from To — 5 sec to To + 40 sec. There is a second double 
peaked structure between To + 95 sec and To + 105 sec. T90 
(15 — 350 keV) is 100 ± 5 sec (estimated error including sys- 
tematics). 

The time-averaged spectrum from To — 14.1 sec to To + 
104.5 sec is best fit by a simple power-law model with a 
power-law index of F = 1.45 ± 0.07. The fluence in the 15 — 
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Figure 1. The mask weighted light curves, at 1 sec time resolu- 
tion, for the BAT energy bands. The double FRED-like structure 
is visible at the start of all of the BAT light curves and a late peak 
is seen at approximately 100 sec. Vertical lines: TgO interval (dash 
lines), T5O interval (dash-dot lines), slew interval (full lines). 

150 keV band is 2.6 zb 0.1 x 10~ 6 erg cm" 2 . The 1 sec peak 
photon flux measured from T — 0.97 sec in the 15 — 150 keV 
band is 1.4 ± 0.1 ph cm - sec -1 . All the quoted errors are 
at the 90% confidence level. 

Using the spectroscopic redshift of this burst z — 3.082 
reported by (jLedoux et al.l H0O6) and a cosmology with 
n m = 0.3, Q A = 0.7, and H = 65, we find E iso (1-1000 keV 
in the rest frame) to be Tlx 10 53 erg. This is based on an ex- 
trapolation of the BAT power-law fit into the corresponding 
observer energy band. 

2.2 XRT Data 

Observations using the Swift X-Ray Telescope (XRT) began 
73.6 sec after the trigg er. It found a bri ght, variable, uncat- 
alogued X-ray source (|Page et al.ll2006|). The XRT position 
enhanced by UVOT astrometrv (|Goad et a"i1l2007bl ) was at 
RA(J2000) = 21:58:50.46, Dec(J2000) = -22:29:47.3, with 
an estimated uncertainty of 1.6" (90% confidence radius). 

The XRT light curve shows three flares peaking at ap- 
proximately 97 sec, 175 sec, and 263 sec after the BAT trig- 
ger. These flares are superposed on a decaying continuum 
with a decay index en = 1.09 ±0.04. At 714tg| sec after the 
BAT trigger the slope flattens, becoming 0:2 = 0.41 ± 0.03. 
There is a second break at 12200^35° sec after which the de- 
cay becomes ay = 3.29tg 1 • The X-ray light curve of GRB 
060607A is shown in Fig. |1 

2.3 UVOT Data 

The Swift Ultraviolet and Optical Telescope (UVOT) ob- 
servations began 63.6 sec after the trigger with a 10 sec 
settling mode exposure with the V filter, followed by finding 
chart exposures of 100 sec with the White (160 — 650 nm) 
filter starting at T p + 73.6 sec and then 400 sec V filter 
jOates et all l200d ) . The settling exposure and each of the 
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Figure 2. The XRT X-ray light curve (0.3-10 keV), WT (black), 
PC (grey). This plot includes all the XRT observations of GRB 
060607A. 

two finding chart exposures are taken in "Image & Event" 
mode, so in addition to the full-frame integrated image, 
event-by-event data is also available at 11.0322 msec maxi- 
mum time resolution. 

A bright afterglow was detected in autonomous ground 
processing software, at RA = 21:58:50.40, Dec = -22:29:46.7 
(J2000) with a la error radius of approximately 0.5". This 
position was 4.7" from the center of the refined XRT er- 
ror circle and is 0.95" from the XRT position enhanced 
by UVOT astrometry. The estimated initial White magni- 
tude in the finding chart image was 15.7 with a la error 
of « 0.5 mag. The optical afterglow was also detected in 
White, V, B, and U filters (Fig[3]). The non-detections in the 
ultraviolet (UV) ban ds are consistent wit h the spectroscopic 
redshift of z = 3.082 l|Ledoux et al.ll2006l ). We have also ex- 
tracted 5 sec binned light curves from event-by-event data 
to investigate the short time variability of the UVOT after- 
glow and to compare it to the XRT and BAT light curves 
(Fig(4]). Throughout this work, all the investigations of the 
UVOT data up to ~ 600 sec is based on this light curve. 

We performed photometry on each UVOT exposure us- 
ing a circular aperture with a radius of 2" centred on the po- 
sition of the optical afterglow. This radius is approximately 
equal to the full-width at half-maximum (FWHM) of the 
UVOT point-spread function (PSF). The PSF varies with fil- 
ter and with the temperature of the telescope, so we did not 
match the extraction aperture to the PSF for each exposure. 
The PSF FWHM, averaged over the temperature variations, 
ranges from l."79 ± 0.05" for the V filter to 2.17" ± 0.03" 
for the UVW2 filter. The background was measured in a sky 
annulus of inner radius 17.5" and width 5" centred on the 
afterglow. 

Aperture corrections were computed for each exposure 
to convert the 2" photometry to the standard aperture radii 
used to define UVOT's photometric zero points (6" for UBV 
and 12" for the ultraviolet filters). Six isolated stars were 
used to compute the aperture correction for each exposure. 
The RMS scatter in the mean aperture correction for a single 
exposure was typically « 0.02 mag. The RMS scatter for 
each exposure was added in quadrature to the statistical 
error in the 2" magnitude to obtain the total la error in 
each point. 

Since the UVOT is a photon-counting device it is only 
able to record one photon per detector cell during each read 
out. This results in coincidence losses at high count rates. 
For very high count rates, corresponding to V < 13.5, these 
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Figure 3. UVOT light curves in White, V, B, and U extend- 
ing out to where only 3-sigma upper limits are seen. They are 
calculated from co-added exposures. The background was mea- 
sured in annulus around the source position as explained in the 
text. The first imaging exposure from ~ 73 sec. to ~ 172 sec in- 
cludes the rising part of the light curve, and therefore the rise 
of the optical flux is only observed in event-by-event data of the 
UVOT (see FigQ. Note also the episode of flattening, or possible 
rcbrightcning, between ~ 1000 sec and ~ 2000 sec. 
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Figure 4. The light curves, at 5 sec time resolution, for the BAT 
(upper), XRT (middle), and UVOT (bottom). The prompt double 
FRED-like structure as well as later spike at ~ 100 sec are visible 
in the BAT light curve. The XRT light curve shows probably 
a minor at the beginning of observations following by the flare 
at ~ 100 sec detected also by the BAT, then one minor flare 
at ~ 170 sec, and one major flare at ~ 260 sec. In the UVOT 
light curve circles correspond to the White filter and bars to V 
filter data points. The UVOT V band count rates are weighted 
by a factor of ~ 3. The optical light curve shows some variations 
at the time of flares, otherwise it is dominated by a continuous 
component peaking at ~ 200 sec. 



losses are significant and can dramatically affect the pho- 
tometry, so coincidence loss corrections must be made. We 
have corrected all of our data for coincidence loss, although 
for most of the observations of the afterglow this correction 
is negligible. 

The values quoted in Table [1] are not corrected for the 
expected Galactic r eddening of E B -v = 0.03 ± 0.02 mag 
ijSchlegel et al.| [l998') . This reddening corresponds to Galac- 
tic extinctions of Au — 0.15, Ab = 0.13, and Ay = 0.10. 
The Galactic extinctions in the UVOT ultraviolet filters 
were calcu lated using the Milky Way extinction law from 
i|Peilll992r i. The ultraviolet extinctions are Auvwi = 0.23, 
Auvm2 — 0.29, and AuvW2 = 0.21 and Awhite = 0.13. 

The afterglow was detected by UVOT from 63.6 sec 
after the BAT trigger until it faded below detectability at 
approximately 2 x 10 4 sec post-trigger. All significant detec- 
tions as well as 3a upper limits are listed in Table Q] and 
shown in Fig. [3] 

There is not a unique way to fit light curves spe- 
cially when they do not have a simple power-law be- 
haviour. For the PROMPT tele s cope o b servations of GRB 
06060 7A l|Nvsewander fc Haislid 120061 ; iNvsewander et all 
2007) authors have used a complex expression including 
absorption terms and power-law rising and falling terms 
around each of the features in the light curve. Here we are 
mainly interested in the origin of the dominant component of 
the optical light curves. Therefore, a simple power-law sepa- 
rately fit on the rising and declining segments of the optical 
light curves is adequate. In Sec[3]we also fit the optical light 
curve by adding a component proportional to the prompt 
gamma-ray to investigate the contribution of the prompt to 
optical emission. 

For the rising section of the optical emission from 
~ 100 sec to ~ 170 sec only White filter observations are 
available. The light curve is fit by a rising power-law with 
a = 2.3 ± 0.3. The relatively large uncertainty is due to the 
varying components during flares in this time interval. The 
decay slopes are listed in Table [5] The weighted mean decay 
index is a = 0.9±0.06 (lcr error) for all times after w 600 sec 
after the BAT trigger. It is consistent with what is seen in 
other optical afterglows before the jet break occurs (Oates, 
et al. in preparation). 

The decay slope is the same before and after the fluc- 
tuation in the light curve between 1000 sec and 2000 sec. 
This f eature was also observ ed by the PROMPT tele- 
scope ([Nvsewander et al1l2007l ) and by the REM telescope 
l|Molinari et al. 2007 ), and is probably the optical emission 
from a weak flare at about 1000 sec after trigger, see Fig. 
[2] Another feature was observed by the PROMPT telescope 
between ~ 3000 and ~ 4000 sec. However, a gap in the Swift 
data from ~ 1500 sec to 5000 sec does not allow to see if it 
is related to an X-ray flare. 

There is also a faint feature in the UVOT light curve 
that at first sight can be considered as noise: A double peak 
at the time of maximum flux. T his feature is independen tly 
observed also by the PROMPT iNvsewander et al.ll2o"o7r ) in 
B and by the REM (jMolinari et al] 120071 ) telescopes in H 
filter, and therefore most probably is real. It can be related 
to the flare at To + 179 sec. 



GRB 060601 A: A GRB with Bright Asynchronous Early X-ray and Optical Emissions 

Table 1. Swift UVOT Photometry of GRB 060607A Table 2. Swift UVOT Decay Slopes of GRB 060607A 
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that a fast varying component is superposed on the domi- 
nant continuous component of the optical emission. In order 
to determine the contribution of this fast varying compo- 
nent to the optical afterglow during early UVOT observa- 
tions, we performed an analysis similar to that of Vestrand 
et al.(2005). We assumed that the UVOT light curve could 
be represented by a component proportional to the gamma- 
ray component in the same time interval: 



Fp(t) = CpF-y(t) 

and a continuous component of the form: 
t-to 



F a (t) = C a 



to 



exp 



t-to 



(1) 



e(t-to) (2) 



where to is an arbitrary initial time, r is the timescale for the 
rise of the optical emission, s is the power-law decay index, 
and C p is the ratio of the UVOT fast varying component in 
White filter (the only filter with simultaneous BAT detec- 
tion) to the BAT 15 — 150 keV flux, and C a is the amplitude 
of the continuous component of the optical emission. The 
step function O is added to restrict this equation to t ^ to. 

Our best fit is shown in Fig. [5] The best fit has 
to = 37.9 ± 2.3 sec, r = 209 ± 5 sec. The power-law index 
has been fixed to s = 2. The ratio of the fast varying 
component to the continuous afterglow in the UVOT data 
is C p /C a = (6.25 ±8.20) x 10~ 2 . This suggests that the fast 
varying optical emission makes a small, if not negligible, 
contribution to the observed optical light at the time 
of the flares. This is in contrast to what was found for 
some of other GRBs with simultaneous ga mma ray, X-ray, 
and optical observations: GRB 041219 A (IVestrand et al.l 



120021) . GRB 050820 A (|Vestrand et all 120051) GR B 



first column shows the instant corresponding to the middle 
of the exposure since trigger time. 

3 ANALYSIS OF TIME VARIABILITY 

The brightness of the early X-ray and optical emissions of 
GRB 060607A together with good time resolution of the 
Swift on-board instruments permit the investigation of the 
relation between various features observed in different en- 
ergy bands. This is very important for identifying the re- 
lated physical processes and their modeling. In this section 
we investigate the correlation between features. 

As mentioned in the previous sections, it is very clear 



060526 (|Dai et all 120071 ). GRB 061121 (|Page et all l2007f l. 
where the prompt/flare component of the optical light 
makes a significant contribution to the total optical light. 
The possible reason can be the fact that for these bursts 
BAT had triggered on a faint precursor that produced a 
faint continuous emission before the occurrence of the main 
gamma-ray peak. In the case of GRB 060607A the main 
peak was at trigger time and had a significant tail emission. 
The time varying component here is due to fainter flares at 
later time, and therefore less significant than the remnant 
of the prompt emission. 

To see whether the X-ray and the optical light curves 
correlate or more precisely if there is any trace of the flares in 
the UVOT light curve, we removed the smooth component of 
the rising section in White filter using the Sa vitzky-Golay 
smoothing algorithm (|Savitzkv fc Golavl Il964l ) , and corre- 
lated the residue with the XRT light curve in 0.3 — 10 keV. 
The feature in the rising part of the White filter light curve 
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Table 3. Spectral Fits to the X-Ray Data 
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Time since BAT trigger (s) 
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Figure 5. The UVOT White and V-band data are fit by the two- 
component model (dashed line) described in Sect. \3\ The data 
from the two filters have been arbitrarily shifted so that they are 
aligned. 



correlates with the flare at To + 98 sec in the XRT with a 
lag of ~ 5 sec. The trace of the flare at To + 260 sec is less 
evident. As mentioned in Sec. I2.3l there is a break in the V- 
band light curve (Fig. 2} with a lag of ~ 5 sec around this 
time, but it is less significant than the first flare. 

As for the common gamma-ray and X-ray flares, Fig. [4] 
shows that they are very close and correlate with each other 
with a lag of < 1 sec with 15 - 300 keV BAT band. There is 
also a weak evidence of a peak close to To + 260 sec and the 
correlation between this section of the BAT and XRT light 
curves gives a lag of ~ 12 ± 5 sec. 



BROAD-BAND SPECTRAL AND LIGHT 
CURVES VARIABILITY 



We fit the X-ray spectra using XSpec-12 (|Arnaudl fl996). 
The model used was a power law spectrum with variable 
hydrogen column density in the host galaxy and a Galactic 
hydrogen column density fixed at N^ IW — 2.67 x 10 20 cm -2 . 
We adopted a redshift of z = 3.082 for this burst and its host 
galaxy. The photon counting data before the break - mainly 
the shallow slope regime - are well-fit by a single power-law 
with T = 1.59±0.06 and an excess N H of 5.7t 2 2 9 7 X 10 21 cm" 2 
at the redshift of the source. After the break the spectral 
slope is r = 1.73 ± 0.08. There is no evidence for a change 
in the column density across the break at ~ 12200 sec. The 
brightness of the X-ray afterglow permits the spectrum to 
be determined before, during, and after flares. From Fig.[5]it 
is clear that strong flares overlap the continuum X-ray emis- 
sion. For this reason we have selected 4 intervals of 10 sec 
during which the flares or the continuum are expected to 
be the dominant contributors to the X-ray lightcurve. The 
duration was selected such that there were enough events 
at later times - smaller fluxes - to obtain a statistically sig- 
nificant spectrum index and Nh- The first and the last in- 
tervals belong to the continuum in the region usually called 
tail emission. The last interval was positioned to be as far 
as possible from the transition region between tail emission 
and the shallow slope regime. The other two intervals are 
positioned near the peaks of the main flares. They allow 
to measure spectral differences between maximally different 
features of the light curve. The results are summarized in 



XRT-Flares 



Time + T (sec) T 



7V H (10 22 cm- 2 ) 



Before 1 st flare 79.5 - 89.5 

Peak 1 st flare 91.5 - 101.5 

Peak 2 nd flare 257.5 - 267.5 

After 2 nd flare 394.5 - 404.5 



+0.26 



2 09 

z.u»_ 23 

1.34 ±0.13 
-,+0.13 
-0.12 
-1+0.42 
-0.45 



1.60Z 

2.00j 



-1.2 
-0.9 

< 0.92 

1+1.0 
-0.8 

< 0.76 



i.r 



l.r 



[^Vh are at the redshift of the source z = 3.082. Galactic Nh is 
2.6 ± 0.05 x 10 20 cm~ 2 . 



Table El We discuss them later in this section along with 
the evolution of the X-ray hardness ratio. 

There are coincident BAT and XRT flares at approxi- 
mately 100 sec, so we performed a joint spectral fit to the 
BAT and XRT data during the flare. The best fit is obtained 
using the same model as described above for the XRT-only 
fits. It has F = 1.63 ± 0.05 and Nh is consistent with the 
value that was determined using only the XRT data. Chang- 
ing the model to a broken power law or a cut-off power 
law does not improve the fit. The joint BAT+XRT spectral 
slope is consistent with what was found for the XRT data, 
so we conclude that the X-ray flare is an extension of the 
flare seen in the gamma-rays. Fig. [6] shows the BAT light 
curve extrapolated to the XRT band using the mean of the 
BAT and pre-flare WT spectra, along with the XRT light 
curve. In addition, it shows power-law fits on the different 
segments of the continuum component of the light curve. 
There is a smooth power-law decaying component, from the 
time at which the first gamma-ray peak begins to fall at 
~ To + 5 sec, until the end of the last major flare in X- 
ray at ~ To + 700 sec. The decay slope of this segment is 
~ 1.17. Flares, in both BAT and XRT, are superposed on 
this continuum. It is evident that flares observed by both 
instruments are in close relation, and therefore X-ray flares 
have the same origin as the prompt gamma-ray emission, 
presumably the internal shocks. As mentioned in Sec[31 the 
correlation of XRT and UVOT light curves shows an ex- 
cess of optical emission during the flares. In conclusion, the 
late flares are simultaneously observed by all instrument on 
board of the Swift. 

A confirmation of the hypothesis that flares have the 
same origin as the prompt gamma-ray can come from a 
direct comparison of their emission properties. There is a 
well known relation between the width of auto-correlation 
of peaks and the energy bands first discovered by Fenimore, 
et al.(1995) - the lower the energy band, the larger the width 
of the auto-correlation of a peak. We tested this relation 
for the BAT peaks and for the flare at ~ To + 98 sec ob- 
served by both BAT and XRT. In all cases peaks have a 
wider auto-correlation for lower energy bands. Another sug- 
gested relation is the proportionality of the waiting time 
during a quiescent period and subsequent burst of radi- 
ation, e vidence for a sort of accumulation of energy be- 
haviour dRamirez-Ruiz fc Merlonl 2001; Ra mirez-Ruiz et al.l 
l200ll ; lNakar fc Piranll2002D . We can not confirm this relation 
for GRB 060607A. In fact, it seems that the amplitude of 
fainter peaks following the main peak at trigger time de- 
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Figure 6. The BAT light curve (black) extrapolated into the 
XRT band using the mean of the BAT and pre-flare WT spectra, 
and the XRT light curve (grey). The red curve shows the fit on 
the continuum component of the lightcurve from To + 5.41 on. For 
To + 5.41 to T + 635.96 the slope is a = 1.17 ± 0.04, consistent 
at lcr level with ai in Sec l2.2l For the other segments the results 
are the same as what is reported in Sec l2.2l 

creases according to a power-law, somehow shallower than 
continuum emission, regardless of the quiescent interval be- 
tween them, see Fig. [6] 

In order to estimate the shortest variability timescale 
of the prompt and flares we use their Fourier transform 
and auto-correlation. The shortest variation timescale dis- 
tinguishable from a white noise for the prompt gamma-ray 
from To — 5 sec to To + 95 sec is ~ 2.5 ± 1 sec in all bands. 
For the flare (peak) at ~ To + 98 sec, the minimum variabil- 
ity distinguishable from noise is close to the main peaks, 
~ 3 ± 1 sec, using the total light curve for the BAT in 
15-300 keV and for the XRT data in 0.3-10 keV. The min- 
imum variability of the flare at ~ To + 260 sec is ~ 7± 5 sec, 
longer than previous peaks. Variation in the emission can be 
due to the inhomogeneity in the shells and/or their limited 
size. In this case one expects that with the expansion of the 
fireball or what rests from it, the density variation dilutes. 
Therefore, what we expect is an increasing minimum time 
variability with time, which is exactly what we are observ- 
ing. However, as with expansion the signal becomes fainter 
and the variations become more difficult to observe, the in- 
crease in minimum variability scales is a consequence of both 
effects. 

As the speed of the passage of the shock front through a 
shell is limited to the speed of light, these timescales can be 
translated to the distance traveled by a relativistic ejecta. 
Therefore, they constrain the initial size and/or variability 
scale of the fireball, and the distance to the central engine to 
> 3 x 10 10 cm in the rest frame of the engine. This is consis- 
tent with the estimation of the model explained in the Ap- 
pendix|X]and applied to the data in Sec(5] The timescales are 
also another confirmation of the same origin for the prompt 
emission and flares seen in X-ray and optical bands. The 
increasing timescale of variations is consistent with the ex- 
pansion of what rests from the prompt shell (fireball) after 
internal shocks and its coalescence with other shells. 

As for the continuum emission, there is no evidence of 
a steep initial decline as would be expected from a high- 
latitude emission. One possibility is that the external shock 
phase with ISM/surrounding material began very early and 
smeared the high latitude emission. Another possibility is 



that the continuous emission and the preceding peaks in the 
prompt emission have the same origin; with the continuous 
emission due to the decaying tail of the emission from the 
prompt shell after its main collision. As explained earlier 
in this section, a power-law fits the extrapolated BAT light 
curve to the XRT energy band. The smoothness of the joint 
lightcurve is a likely evidence that the initial smooth, some- 
times very steep decay observed in the X-ray light curve of 
many GRBs is directly related to the prompt emission and 
hence the term tail emission is a correct expression for this 
regime. As for the physical processes involved, it is possible 
that the magnetic field in the coalesced shells has a relatively 
long lifetime and electrons are accelerated and support a 
synchrotron emission at a lower rate well after the end of 
the collision between shells. See also Sec. [5] for more details. 
There is however a caveat in this argument. The extrapo- 
lation of the BAT light curve to the XRT uses the average 
spectrum slope observed by these instruments in their cor- 
responding bands. Therefore, the smoothness of the joint 
light curve may be due to the way it is calculated, thus the 
argument about common origin of the prompt and tail emis- 
sion becomes doubtful. Nonetheless, we show in Appendix 
[B]that if the time evolution of the BAT spectrum index is 
taken into account, it is highly improbable to obtain such a 
smooth common light curve due to averaging or by chance, 
and one should observe a deviation between extrapolated 
BAT light curve and what is observed by the XRT if there 
is not an intrinsic relation between the prompt and the tail 
emissions. 

Fig. [7] shows the evolution of the X-ray hardness ra- 
tio: HR X = (7(1.5 - 10 keV)/C(0.3 - 1.5 keV), where C 
is the count rate uncorrected for the absorption. Both of 
the X-ray flares are significantly harder than the underlying 
decay, consistent with the spectral indices reported in Ta- 
ble [3] Moreover, comparing the spectral index of the spec- 
trum of the first and the last time intervals reported in this 
table, and their hardness ration shown in Fig. [7] it seems 
there is no evidence for large spectral evolution in the con- 
tinuum emission until entering to the shallow slope regime. 
The hardness ratio during each flare tracks the luminosity; 
that is, the hardness increases as the flare brightens, and 
decreases as the flare fades. Another interesting observation 
from this plot is a gradual rise of the hardness ratio from 
the beginning of the shallow regime to a roughly constant 
plateau that does not decline even after the sharp break at 
~ 1.2 x 10 4 sec. If the break was achromatic such a be- 
haviour was expected, however the break is not achromatic. 
It can be due to a decrease in the number of emitters - accel- 
erated electrons - and hardening of their spectrum, and/or 
a stronger magnetic field. In Sec. [S] we discuss these issues 
more in details. 

The fluence in gamma rays (15-150 keV) during the first 
pulse (T - 24 to +12 sec) is (1.38 ± 0.05) x 10" 6 erg cm" 2 
while the corresponding fluence between 24 sec and 102 sec 
is (1.22 ± 0.05) x 10~ 6 erg cm" 2 . There is no flux in the 
15 — 150 keV band after 102 sec. Therefore, unlike the first 
flare the second XRT flare has no (or a marginal) BAT coun- 
terpart. This suggests that the two XRT flares are produced 
by distinct shells colliding with the remnant of the prompt 
shell. 

It seems that both main flares were preceded by fainter 
but harder flares, see Fig[6] This repetition can be explained 
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Figure 7. The upper and middle panels show the count rate mea- 
sured by the XRT in 1.5 — 10 keV and 0.3 — 1.5 keV respectively. 
The lower panel shows the evolution of the hardness ratio, the 
ratio of the flux measured in the 1.5 — 10 keV band to that mea- 
sured in the 0.3 — 1.5 keV band. Darker points are observations in 
WT (Window Timing) mode and lighters in PC (Pulse Counting) 
mode. 



if two late shells successively pass through the remnant of 
the prompt shell. In this case the observed emission is the 
tomography of the prompt. The small and large successive 
flares can be the result of collision between the late shell, 
slower shocked material, and a leftover unshocked mate- 
rial from the prompt shock with a larger Lorentz factor 
that moves ahead of the first component. A similar con- 
figuration of fl ares is also visible in t he X-ray light curve 
of GRB 070107 l|Stamatikos et aDl2007h . Assuming that the 
late shells have roughly the same dynamical properties - 
Lorentz factor, density, etc. - the softness of the second 
flare could be due to the expansion and slowdown of the 
unshocked remnant. The expansion of the late shell can also 
be in part responsible for the softer radiation during colli- 
sion, which in this interpretation had necessarily happened 
at larger distances from the central engine. 

The observed X-ray decay slope after flares and up to 
~ To + 600 sec when the shallow regime began, is ax = 
1.09 ± 0.04. The spectral index of the X-rays in the same 
interval is fix = T — 1 = 0.64 ± 0.07. If the emission in this 
interval is due to an external shock with ISM/circumburst 
material, the closure relationships can be used to determine 
the density distribution. We consider two cases: a constant 
density circumburst medium and a wind-stratified circum- 
burst medium. The closure values are given in Table [4] The 
case that gives a closure closest to zero is one with a con- 
stant density medium and the cooling frequency above the 
X-ray band at 600 sec. Further evidence for the cooling break 
being above the X-ray band in this time interval is that 
the optical decay has a slope of a opt = 0.9 ± 0.06, close to 
the X-ray slope at la and consistent with it at 2a level. If 
this scenario is correct, the electron index is predicted to be 
q — p + 1 ~ 3.5, consiste nt with what is seen in many other 
GRBs (|Shen et al.ll2006h . The optical decay does not change 
between ~ 600 sec and ~ 20000 sec after the BAT trigger, so 
the cooling break must be above optical frequencies during 
this period. 

As for what we can learn about the gas and dust content 



of the host galaxy of GRB 060607A, the fitted neutral hydro- 
gen column density in the host along the line of sight to the 
burst is consistent with Nh = 1.6 x 10 22 cm" 2 , which implies 
a hig h extinction in the host galaxy. l|Predehl fc Schmittl 



1995) find N H = (1-79 x 10 21 )A V for the conversion be- 
tween hydrogen column density and extinction in the Milky 
Way. If this relationship holds for the host galaxy, then 
Av = 8.9 mag in the host. While the optical data does 
not rule out such a high extinction there is no evidence for 
it either. If we assume an Nh /Av ratio like that in the SMC 
of Nh = ( 15.4 x 10 21 )ylv (from equation (2) and Table 2 of 
l|Peilll992l ) then A v in the host is 1.0 mag. 

Gas-to-dust ratios similar to that of the SMC have 
been obse rved for several GR B host galaxies, such as GRB 
00030 1C l| Jensen et all ,20011) . GRB 00926 l|Fvnbo et all 
|2001| ). GRB 020124 dHiorth et all ,2003,1. and XRF 
050416 A jHolland et alj|2007h . see also (|Schlegel et al.lll998l ; 
iKahn et al.l ,20061). Detailed investigation of 7 Swiftbursts 
( Schadv et al. 20071 ) shows that only in one case the extinc- 
tion is best modelled by Milky Way gas-to-dust ratio, and 
for other cases the extinction is more similar to SMC. Giv- 
ing the fact that GRB060607A was observed in all the rest 
frame bands redder than Lya suggests that it may be rea- 
sonable to assume that the host galaxy of GRB 060607A 
has a high gas-to-dust ratio. This could indicate that star 
formation in the host is fairly recent and there has not been 
enough time for large amounts of gas to be processed into 
dust. Alternately, a high ratio could be indicative of dust de- 
struction in the vicinity o f the progenitor by t he burst itself 
llWaxman fc Draind |2000| ; iPerna et al.l ,20031 ; IWatson et al.l 
,20071) . 

There is no evidence for a jet break in the optical light 
curve out to ~ 20000 sec after the BAT trigger. This, and the 
isotropic energy of the burst, can be used to put a lower limit 
on the opening an gle of the .jet, and the total gamma-ray en- 
ergy of the burst (|Rhoadslll999l ; ISari et al.|["l999l ; iFrail et all 
2001). Assuming a single jet, the lower limit of the jet open- 
ing angle for GRB 060707A is: 



> 0.161 



z+1 



3/8 frnj^ 1 ' 



1/8 



rad, 



(3) 



where is the efficiency of converting energy in the 
ejecta into gamma rays, and n is the particle density in 
cm -3 . The corresponding energy in gamma rays, using Ei so 
obtained in Sec. 12.11 in the rest frame of the engine and 
corrected for the beaming is E 1 ^ 3.3 x 10 49 erg. 



5 MODELLING 

In this section we apply a reformulation of the inter- 
nal/external shock model for the prompt and afterglow to 
the data in order to estimate some of the parameters of 
the burst. This formulation is based on a simplified ultra- 
relativistic radiative shock model with one synchrotron emit- 
ting shocked layer but more detailed parametrization of the 
physical processes and their time variation (Ziaeepour, in 
preparation). A summary of the model and its main results 
are given in the Appendix IA1 
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Table 4. This Tabic lists the closure relationships for various as- 
sumptions about the location of the cooling frequency and the 
nature of the circumburst medium. A closure value of zero indi- 
cates agreement with the predictions of each case. 



Model 



Environment Closure 



Value 



ISM 
Wind 

ISM 
Wind 



■3/2/9 

■3/2/3 + 1/2 
•3/2/3- 1/2 
■ 3/2/3 - 1/2 



-0.13 
+0.63 
-0.37 
-0.37 



5.1 Methodology 

A major difficulty in understanding the behaviour of GRBs 
is that the main ingredients of the shock i.e. the elec- 
tric and magnetic fields, and the distribution of electrons 
varies in a complex manner with time during the evolu- 
tion of the microphysics and the dynamics of the shock 
llWiersma fc Achterberd [iooj; iBednarz fc Ostrowskil 1 199rj; 
Waxman fc Draind l2000t IReville et all 120061 ; IRieger et ail 
20061 ). The commonly used power-law parametrization with 
constant coefficients and indexes is not able to explain 
the complex behaviour of quantities and thereby the syn- 
chrotron emission. Adding the microphysics of the shock to 
the formulation of the shock dynamic is also too difficult 
as in most cases there is no analytical expression for their 
evolution, or for the evolution of the electric and magnetic 
fields or distribution of electrons. As we will explain in more 
details, we found that the best way to estimate parameters 
and to explain the behaviour of the light curves is to divide 
them into sub-regimes that can be explained separately by a 
simple power-law parametrization. On the other hand, this 
simplified model still requires 12 parameters and thus fit- 
ting such a complex model to the data is not trivial, even 
after sub-dividing it into separate regimes. The degeneracy 
between parameters and the instability of numerical fitting 
can lead to confusing results. Moreover, at the current level 
of our knowledge about GRBs and the precision of avail- 
able data, even a rough estimation of parameters permit a 
better understanding of the nature of these elusive objects 
and can be considered as an achievement. Therefore, rather 
than fitting the data, we determine the predictions of the 
model through numerical simulations with prefixed param- 
eters. Then, we estimate the set(s) of parameters that best 
reproduce the behaviour of the data. 



5.2 Modelling and interpretation of the BAT data 

To interpret the light curves of the GRB 060607A, we be- 
gin with the hardness ratios of the BAT bands. The hard- 
ness ratios of the first main peak are shown in Figs. [S] The 
initial rapid rise is consistent with an exponential grow- 
ing to a saturated value of e e and es, respectively, of the 
fraction of kinetic energy transferred to accelerated elec- 
trons and to magnetic field. This is also in acco rdance with 
simu l ations of the formation of mag netic field l|Yang et al.l 
1 1994 Wicrsm a'fc Achterberd I2004T ). parallel instabilities 
I Reville et al.l 2006) in the relativisti c shocks, and acceler- 
ation of particles in a Fermi process (jBednarz fc Ostrowskil 




time - 1.7135e+08 (sec) 

Figure 8. Hardness ratios of the main peak. From top to bot- 
tom: HR 10 ee C*(25 - 50 feeV)/C(15 - 25 keV), HR 2 i = C(50 - 
100 fceV)/C*(25-50 keV), and HR 32 = C(100-300 keV)/C(50- 
100 keV). For calculation of the ratios data from the 64 msec 
binned light curve are used and rebinned by a factor of 30. 



Il99rj ; IRieger et all 12006^1. After this transient time, fields 
settle into a quasi-stationary regime in which they decline 
according to a power-law. Comparison with simulations 
show that for this burst r\ = 2a e + as/2 ~ 3.5 ± 0.5 for a 
power-law electron number distribution index q ~ 3.5 ± 0.5, 
where a e and as are time/radius dependence of e e and es, 
respectively. The value of q ~ 3.5 (or equivalently p ~ 2.5) is 
close to p ~ 2.2 — 2.3, the universal value suggeste d by sim- 
ulations (|Kirk et al.ll2000l : lAchterberg et al[|200ll '). On the 
other hand, the distribution of p obtained from the observa- 
tions of GRBs and other relativistic sources does not show a 
universal averag e and has a relativ ely large deviation around 
the mean value (|Shen et al.ll2006h . This can be the evidence 



3 For simplicity in all the simulations discussed here we use a 
power-law with constant index distribution for electrons. 
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Table 5. BAT lags 



Band 


Lag (ms), entire burst 


0-2 


621.90 ± 57.0 


1-2 


252.18 ± 25.7 


3-2 


-556.73 ± 57.0 



for more complex outflow behaviour such as: a complex ge- 
ometry of the magnetic field and shock, plasma currents and 
turbulence, non-Fermi acceleration mechanisms, or domi- 
nance of the Poynting flux. In conclusion, although the p 
value obtained here is typical, it can not give more informa- 
tion about the underlaying processes. The estimated value 
of characteristic synchrotron energy of the electrons with 
minimum Lorentz factor V is £ m = hv m ~ 12(1 + z) keV 
~ 50 keV in the rest frame of the engine. 

These estimations are based on the hardness ratios 
HRu) and HR21 ■ The highest energy band hardness HR32 = 
C(100-300fceU)/C(50-100fceU) is much noisier than H R w 
and HR21. However, it seems that the electron spectrum 
should be much steeper to explain the observed low value 
of HR32 ~ 0.4 in place of the expected value of ~ 0.65 
from simulations for q and r\ as mentioned above. We in- 
terpret this as an upper limit cutoff or steepening in the 
electron di stribution. Simulations o f electron acceleration 
by shocks (|Amano fc Hos hino 2007) also confirm a steeper 
distribution at high energies. A change in the spectrum 
of electrons has been also invoked as an explanation for 
the chromatic break in the X-ra y and optical afterglows 
l|Wei fc Lull2002l : lMisra et al.ll2007h . The observation of such 
a break/steepening in electron energy distribution during 
the prompt emission, when the shock is much stronger, 
shows that the popular assumption of a power-law distri- 
bution is too simplistic and far from reality. The estimated 
range of parameters includes both uncertainties and degen- 
eracy between them. The initial (maximum) value of the 
hardness ratios depends mainly on q, their decreasing rate 
depends on 77 and coefficient F defined in IA3I and there is 
a significant degeneracy between these quantities. 

From the slope of the hardness ratios in the quasi-steady 
regime and the lag between the BAT energy bands, we can 
estimate the coefficient F and the distance to the central 
engine. The lags between the BAT bands are presented in 
Table [5] Considering the peak time after the initial rise as 
the beginning of the power-law decline regime for e e and 
e_B, and by comparing lags with simulations, we conclude 
that F ~ 12 - 18, Ae 2 ~ 6 x 10~ 6 , and Ae 12 ~ 3 x 10" 6 , 
where Aejj is the difference between e defined in Eq. IA3I at 
peak time for bands i and j. Including uncertainties both 
in data and comparison with simulations, we estimate the 
initial distance to the central engine as ro ~ 10 12 cm. This 
is roughly in the mi ddle of the distance ran ge predicted for 
the internal shocks l|Rees fc Meszarosl[l9 94n. 

Knowing F and ro, we can also estimate the relative 
Lorentz factor of the shells. In the radiative shock model 

4 In all the calculations here we assume that the bulk Lorentz 
factor of the fireball -fbulk 3> 1 and fibulk ~ !• 



assumed here colliding shells coalesce. Therefore, at the end 
of the shock, 7 — * 1. Using Eq. IA3I 70 ~ 1 + Ftf where 
€f corresponds to the value of e at the end of the coales- 
cence. It is not very evident what time should be used to 
determine e/, because substructures/overlapping peaks can 
be due to separate shells or from density inhomogeneities in 
the same shell. Moreover, the limited sensitivity of the BAT 
can smear the real collision time. Nonetheless, if we assume 
that the detected duration of the peaks corresponds to the 
main part of the collision, an estimation of 70 ~ 1.5 — 3.5 
for overlapping peaks from ~ To — 5 sec to ~ To +40 sec can 
be made. This includes uncertainties in the parameters and 
the duration of the collision. The bulk Lorentz factor of the 
ejecta however can not be determined from the prompt with- 
out knowing all the parameters such as the density of shells, 
magnetic fields, etc. The reason is the fact that the physics 
of shock depends only on the relative Lorentz factor, and the 
detected radiation is just boosted by the bulk Lorentz factor 
to the observed energies. Other BAT peaks are too weak to 
permit a detailed analysis of their corresponding shock. 

5.3 Modelling and interpretation of the XRT data 

If we neglect flares that are superimposed on the early XRT 
light curve shown in Fig. [2] the power-law component from 
~ To + 73 sec to ~ To + 700 sec has a slope of ai ~ 1.09 and 
spectrum index of /3i ~ 0.64. These indices do not satisfy the 
relation qi = /3i+2 for high latitude emission. Therefore, we 
interpret this section of the X-ray light curve as the tail emis- 
sion from energy dissipation in the prompt shell after the 
termination of the coalescence. Investigation of the hardness 
ratio of the X-ray bands C(1.5 - 10 fceV)/C(0.3 - 1.5 keV) 
shows that at this time the ejecta has significantly cooled, 
consistent with the results of Sec. [4] In fact, we can go fur- 
ther and determine £ m at the beginning of the observations 
after Swift slewed to this burst. 

Assuming synchrotron radiation as the source of the 
observed X-ray, according to our model, £ m evolves ao 

^)=Mro)(£f g)* (4) 

After the coalescence of the two shells 7 = 1, and there- 
fore at ~ To + 73/ (z + 1) sec in the source frame when the 
coalescence of the shells has been already finished, the sec- 
ond term in Eq. |?]is fixed to (l/7o) 5 ^ 4 - Using ro, £ m (ro), 
and 70 ~ 2.5 - the mean value in the range derived from 
BAT data - we find that £ m has been reduced by a factor 
of ~ 0.026 to £ m ~ 0.3(1 + z) keV ~ 1.2 keV . Our simula- 
tions show that the observed X-ray hardness ratio (Fig. [TJl 
of ~ 0.5 ±0.2 just before the first flare is consistent with the 
reduced value of £ m mentioned above, only if the spectrum 
of the electrons also has steepened from q ~ 3.5 ± 0.5 to 
q ~ 5. The steepness of this slope may signify the failure of 
the simplified model without an upper limit or high energy 
break in the energy distribution of accelerated electrons. As 
mentioned above, this assumption is not realistic and can 
compromise the interpretation of the data. 

5 Here we neglect the change in the bulk Lorentz factor. This 
happens mainly during the coalescence. 
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The value we obtain for the decay index of the fields 
r\ ~ 4 is also slightly higher than its value after the steady 
state regime during the main peak in the prompt gamma- 
ray emission r\ ~ 3.5. Both the increase in r\ and electron 
distribution index q are consistent with the cooling of the 
shell in this regime. What is remarkable is that we do not see 
a faster decay of the fields after the collision and coalescence 
of the shells both terminated. This means that the coherent 
bulk and field structures formed during the collision have 
relatively long lifetimes and do not disappear immediately 
after the end of collision. 

Fig. [7] shows the hardness ratio of the X-ray bands 
for multiple flares detected by the XRT. The flare at ~ 
To + 98 sec was also observed by the BAT, and there is likely 
a small but delayed trace of it in the UVOT light curve, as 
explained in Sec. [3] A noticeable difference between the flare 
hardness ratio in the X-ray and the gamma-ray emission is 
that the hardness ratio of the X-ray flares follows, sometimes 
with lag or lead, the X-ray light curve, whereas, the hardness 
ratio of the gamma-ray emission, rises to a plateau and then 
declines slowly until the end of the spikfl A strong correla- 
tion is usuall y observed in flare s where flares that brighten 
also harden ( Goad et al.ll2007ah . This is interpreted as the 
consequence of the spectrum evolution and a decreasing 
break energy. However, in the case of GRB 060607a, a high 
energy break was not observed in the XRT or BAT spec- 
trum. Even without a break, the spectral evolution during a 
flare could still produce this correlation. Nonetheless, rapid 
variation of the hardness ratio of flares in GRB 060607a 
does not seem to be a typical behavi our of all bursts, see for 
instance the case for GRB 061121 (|Page et all 120071 '). One 
explanation could be a strong absorption of the soft band 
photons which flatten the soft X-ray light curve. However, 
the flux hardness ratio curve in which corrections are made 
for the absorption has a very similar behaviour, and there- 
fore absorption can not be the cause. Another possibility is 
that the magnetic field and electron acceleration during the 
collision of the late shells with the remnant of the prompt 
shell did not achieve to rise to the steady state regime and 
decayed exponentially once the late shell either coalesced 
with the prompt or passed through it. Our simulations also 
confirm this possibility. 

As for the physical reasons for such behaviour, one im- 
portant factor can be the he ating of the prompt shell after 
its collision. It is well kn own (Wi ersma fc Achterberdl2004l ; 
IWaxman fe~Dr ainc 2000) that the development of a coher- 
ent magnetic and electron acceleration is weaker when the 
shock medium is hot and particles have a significant momen- 
tum in the direction perpendicular to the direction of the 
boost /bulk movement. As the hardness ratio in the flares 
does not seem to arrive to a steady state regime, and we do 
not know the size and other characteristics of the late shells, 
we can not estimate shock parameters for them as we did 



6 Note that we compare the behaviour of the X-ray hardness ratio 
of the main BAT peak with XRT flares. Although the first XRT 
flare of GRB 060607a was also observed by BAT, it is too weak to 
permit the calculation of a meaningful hardness. Nonetheless, the 
similarity of the profiles in Fig.[8]suggests that the XRT hardness 
which is calculated for energy bands just less than one order of 
magnitude lower should most probably be a continuation of the 
same type of behaviour. 



for the prompt gamma-ray peak. Nonetheless, from our sim- 
ulations using q and r\ similar to the prompt and continuous 
component of the X-ray light curve, we conclude that the 
shock is soft, consistent with weak or non-observation in the 
BAT bands. 

Like other bursts the most difficult part of the X-ray 
light curve to explain is the shallow slope regime. Vari- 
ous proces ses are suggested, such as: refreshed shocks b y 
late shells (IRees fc Meszaro A 1 19981 ; ISari fc MeszaTol l200Ch ■ 
continuous energy injection to the fireball (prompt shel l) 
jDai fc Lulll998l : IRees fc Meszarodliootil ; lYang et al.lll994l ). 
and variation of the microphy sics of what remains from the 
prompt shell l|loka et al.l|2006i ). The first two processes seem 
to have an energy problem; they need mu ch larger efficienc y 
for gamma-ray emission than expected (lloka et al.l I2006T ) . 
In addition, the uniform slope of this regime in all bursts 
does not look like a phenomenon relaying on the random 
injection of late/slow shells. The abrupt break at the end 
of this regime - especially in the case of GRB 060607A - 
needs a sudden stop of both energy injection and radiation 
that seems unphysical. By contrast, a change in the micro- 
physics and energy dissipation of the prompt remnant and 
later shells seems a more reasonable cause. In fact, we see 
such a late time slow evolution of the intensities in the simu- 
lations. However, we do not see an abrupt break. Assuming 
that the late break in the X-ray light curve is not the jet 
break as it is not achromatic, its absence in simulations may 
be due to the limited precision of the numerical simulations 
and the simplicity of the model. Moreover, the fact that 
our analytical approximations are valid only for e < 1, and 
therefore, simulations are limited to this linear regime can 
be another factor influencing the lack of a break in the sim- 
ulations. If this interpretation is correct, parameters of the 
shell are similar to the tail emission with a steeper slope. 
The X-ray hardness ratio is however higher than the tail 
emission see Fig. [7] Energy injection by flares to the prompt 
shell can be the reason for the hardening of the synchrotron 
emission. Another explanation can be the onset of a for- 
ward external shock by the ISM or circumburst material at 
a distance ~ 10 14 cm from the central engine. 

The question which arises here is: At these late times, 
what keeps the coherent magnetic field and the electron ac- 
celeration in the shell going? Do we need a continuous fall 
of circumburst material or arrival of late shells (i.e. energy 
injection) to the shell to keep the radiation level high? Are 
the stock of electrons and residual magnetic field enough to 
keep the low flux of the late time radiation? In the first case 
either the circumburst material should exist at all distances 
after the prompt collision, or there must be a discontinuity 
in the emission between the end of the internal shock and 
the onset of the external forward shock. The same type of 
argument is relevant for the continuous energy injection. No 
discontinuity has been observed by Swift. Only a detailed 
knowledge of the origin of the ejecta, the shock, the state of 
the matter in the remnant after the prompt shock termina- 
tion, and the surrounding material can clarify these issues. 



5.4 Modelling and interpretation of the UVOT 
data 

Finally we try to explain the optical light curve of GRB 
060607A which has interesting behaviour rarely seen in other 
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Figure 9. Estimation of the time evolution of £ m in observer 
frame, from the end of the gamma-ray peak until the optical peak 
according to the model explained in the Appendix lAl 



bursts. In general, optical light curves of Swift bursts within 
the first few hundred seconds after the trigger, have a broad 
range of temporal indices, mostly consistent with decaying 
behaviour. After approximately 500s, all light curves are de- 
caying (Oates et al. in prep). In few bursts including GRB 
06060 7A (see Figs. H and GRB 060418 ilMolinari et all 
120071 ). and GRB 070616 jStarling et al l l2007bl ) an inT 
tial rise in the optical light curves is detected in all fil- 
ter from blue to infrared. A number of interpretations for 
this beha viour have been put for ward: The onset of exter- 
nal shock(|N vscwand er et al.ll2007h . forward shock with ISM 
(Molinari ct al. 2007), deceleration of the forward shock jSarJ 
ll997ft fOates~ et al. in prep). Here we argue that the most 
plausible reason for this behaviour is the gradual cooling 
and energy dissipation in the prompt shell and the entrance 
of Em to optical bands. The relatively shallow slope of the 
early X-ray light curve, a = 1.09 may be the consequence of 
a slow evolution of £ m which permitted the detection of its 
transition through the optical bands. 

Using Eq. [4] and the argument about the evolution of 
7 2 /7o> we can estimate £ m at the peak time in optical light 
curve t pea k ~ To + 200/ (z + 1) sec in the source rest frame, 
and we obtain £ m ~ 4(1 + 2) eV = 16 eV. This corresponds 
to ~ 440 nm, with in the optical blue band for the observer. 
Therefore, we interpret the maximum of the light curve as 
the time where the characteristic synchrotron emission for 
the least energetic electrons enters to the optical bands. The 
approximate time evolution of £ m according to the model 
used here is summarized in Fig[9] 

The rising slope of the optical light curve is a r op B t ~ 
2.3 ± 0.3. According to the model, the slope in this regime 
is q/2 — 1/6. With the estimated q ~ 5 from the X-ray data, 
the observed value is consistent with the interpretation of 
the peak as the passage of £ m to optical bands. Therefore, 
another component - external forward or reverse shock - is 
not required to explain the observations. Although there are 
UV observations for this burst, they do not begin until hun- 
dreds of seconds after the rise. Therefore it is not possible to 
see if they peak earlier than longer wavelengths, as would be 
expected from this interpretation. On the other hand, obser - 
vations by the PROMPT telescope (|Nvsewander et al.ll2007l ) 
in B, g' , r' , and i! are consistent with this interpretation. 
In fact the energy difference between B and i' bands with 
respect to the energy difference between X-ray and optical 
is very small, and therefore the expected time shift between 



their arrival to maximum is very small and difficult to detect 
with the time resolution of ground telescopes. Only UVOT 
in event mode could be able to directly detect the shift, 
but it can not simultaneously observe with multiple filters. 
It may be possible to det ect a small shift in the m aximum 
by fitting the peak, but (|Nvsewander et all 120071) did not 
attem pt this. Nonetheless, in Fig. 1 of ( Nvsewander et all 
l2007t) w hich, includes the H b and data from the REM obser- 
vations dMolinari et al-lfeoOTT ). shows a small delay between 
the p eak time in the H band at 180 ± 6 sec (|Molinari et al.l 
2007) with respect to higher energy bands. We should also 
mention that flares observed in gamma- and X-ray bands 
are marginally distinguishable in the UVOT observations. 
This shows that as expected, the synchrotron emission from 
the flares is much harder than the cooling tail emission, 
and consequently no significant optical excess is observed 
during the X-ray flares. The lack of evidence for a signifi- 
cant external shock in the early low energy bands emission 
makes the claim of d etermination of the bulk Lorentz factor 
jMolinari et alj|2007l ) unlikely, as it is based on the assump- 
tion that this peak is the result of a shock on the ISM. In this 
case jbulh is the same as collision 7 because ISM is roughly 
at rest with respect to the source. But according to our re- 
sults the peak is related to the prompt shock and so the 
ISM/circumburst material had no significant contribution 
in its formation and properties. 

The rest of the optical light curve follows a power-law 
decay at least until ~ To + 25000 sec without any change 
in the slope. The magnitude limits at later times are just 
brighter than the expected fluxes from a power-law decay 
and we can not claim that any late break in optical/UV 
bands has been observed. The flattening or slight bright- 
ening from ~ To + 1000 sec to T + 2000 sec, (see Fig. 
is probably related to the weak X-ray flare that occurred 
during this period. 

Although the number of data points in the light curve 
in this regime is very limited, one can speculate that the 
slope of the decline is slightly steeper in longer wavelengths. 
This can be interpreted as the hardening of the emission, 
probably due to fireball encountering the ISM/circumburst 
material. In other words, it seems that a weak external shock 
starts a few thousands of seconds after the trigger and re- 
freshes the prompt shell. Regarding just optical data, this 
claim seems very speculative. However, the hardness ratio of 
X-ray bands is definitely increasing (Fig[7|. This increases 
the hardness ratio to ~ 1, i.e. £ m moves to the X-ray band 
but with a smaller density for high energy electrons. As the 
shell is already decelerated, a few thousands of seconds later 
at ~ To + 1.2 X 10 4 sec the additional mass of the coalesced 
material from ISM breaks the coherence of the magnetic field 
and the synchrotron emission at higher energy bands begins 
a steep power-law decay due to a break in the spectrum of 
the electrons. This does not affect the emission in the optical 
bands from slower electrons. In this sense one can say that 
the X-ray and optical emissions come from separate com- 
ponents. It is also possible that at least part of the optical 
emission comes from the heated-shocked material without 
the presence of a coherent magnetic field. In conclusion, the 
break is not likely to be associated with the jet break, i.e. 
when 1/fbulk > & jet, and the relativistic collimation is re- 
moved. 
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5.5 Alternative interpretations 

To end the discussion, we consider some of the alternative 
interpretations of the optical light curves. First we consider 
the possibility that the optical emission is produced by a 
forward or reverse external shock, and thus has a different 
origin than the prompt gamma ray and the early X-ray emis- 
sion. Considering first the forward shock case, the peak of 
the optical emission marks the onset of the afterglow where 
the peak time of the optical emission, t pea k, is associated 
with the deceleration time of the outflow tdec, and with the 
optical emission at t > idee, most likely dominated by the 
forward shock with £ op t > £ m . This conclusion is similar 
to the results we obtained by considering the early X-ray 
and optical emissions from the tail emission in the prompt 
shell without evoking a second component and a separate 
origin for the optical and X-ray emissions. In fact, using the 
observed power-law spectrum of the early X-ray, we can es- 
timate the expected optical flux by extending the spectrum 
to the optical bands. What we find is summarized in Table[6] 
The observed magnitudes in this table are corrected for the 
Milky Way absorption which is only ~ 0.1 magnitude. They 
are not corrected for the host absorption because we do not 
have any information about the dust to gas ratio in the host 
galaxy. As discussed in Sec(4] the host of the most GRBs 
have values similar to SMC, Ay ~ 1 magnitude. Although 
the difference between the expected and the observed may 
partially be due to the absorption in the host, we also no- 
tice that discrepancies are time dependent, see Tabl^S] The 
value of T4bs — V ple d varies from ~ 1.5 mag up to ~ 5 mag 
depending on the time interval. The observational errors of 
the observed magnitudes are at most ±0.25 at la or ±0.5 
at 2a. They can not explain the large variations as well as 
differences between expectations and observations. There- 
fore, we conclude that at least part of the difference between 
observations and expectations is intrinsic to the fireball. If 
the early X-ray and optical radiations have different origin, 
then the discrepancies between expectations and observa- 
tions become even larger. This is because we must add the 
contribution in the optical flux of the other component to 
the contribution of the component producing the X-ray. 

We now consider the possibility of a reverse shock as 
the source of the early optical emission. The rising slope 
of the optical light curve at t < t pea k is ~ 2.3, which - 
in the context of the standard reverse shock model that 
predicts a slope of ~ 1/3 - is too steep for the passage 
of £ m through £ op t- The rising slope is a probe of the 
strength of the reverse shock and the observed slope suggests 
that the reverse shock is weak - the dimens ionless parame- 
ter C = ([3£y(47rnm p c 2 )] 1/3 /Ar ) 1/2 r- 4/3 ||Nakar fc Piranl 
12004 ) presenting the strength of the reverse shock is large 
> 2. This lower limit on £ is valid both if the emission at 
t < t pea k is dominated by the forward external shock or if 
it is from a non-relativistic or at m ost mildly-relativistic re- 
verse shock l|Nakar fc Piranll2004D . Therefore, it seems that 
for this burst the contributions of these components are 
small. Nonetheless, if the early optical emission is due to 
a reverse shock, the temporal slope of the rising phase also 
probes the external density profile. A wind-like density pro- 
file would produce a significantly slower rise than a uni- 
form external density. The steep rise seems to favor a uni- 
form external density over a wind-like density profile below 
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and around the deceleration radius rd ec - More specifically, 
a wind-like external density profile could not produce such 
a sharp rise either by a forward or through a reverse shock 
even if the reverse shock is extremely weak (( > 1) and 
the outflow acts as a perfect piston as far as the external 
medium is concerned. 

Another possible explanation is that the optical emis- 
sion is mainly produced by the combination of reverse shock 
and energy injection. The latter is one of the possible causes 
for the slow decay of the optical light curve after the peak 
~ i -1 ' 3 compared to the ~ t~ 2 expected for an idealized 
reverse shock where the outflow has a sharp lower cutoff in 
7. The observation of X-ray flares for this burst is also con- 
sistent with the additional outflow after the prompt shell. 

There is a mino r break in the optica l light (seen more 
clearly in the data of (|Molinari et al.ll2007l ) at ~ To +103 sec, 
roughly at the end of the bright X-ray flares and before 
the small X-ray /optical flare. This may be the evidence of 
the break in the reverse shock after a relatively flat regime 
following the initial rise l|Zhang fc Meszarosf2 002 ) . However, 
for this argument to be viable, the reverse shock should be 
significantly brighter than the forward shock. Regarding the 
arguments above about the weakness of the reverse shock, it 
does not seem that we can associate this break to the reverse 
shock. 

We should also mention that the above break can be 
also explained as a hardening of the radiation at the arrival 
of a new shell responsible for the weak flare in XRT data (see 
Fig. [2j at ~ To ± 1260 sec. A similar break is also visible in 
the UVOT light curve during the flare at ~ T ± 260 sec, 
followed by a slightly steeper fall of the optical light curve 
most probably due to energy injection and hardening of the 
emission (see Fig. 2}. Therefore, considering all the argu- 
ments together, the evidence for a dominant reverse shock 
is weak. 



6 CONCLUSION 

After years of investigation and especially with the multi- 
band data from Swift of more than 200 bursts, it is clear 
that complex physical processes are involved not only in the 
physics of the progenitor but also in the collisions between 
ejected material and the surroundings that at first sight seem 
simple. 

In this work we applied the internal/external shock 
model as the source of observed emission to GRB 060607A, 
one of the most peculiar gamma-ray bursts observed by 
Swift. We tried to understand and interpret the acquired 
data and related processes partly qualitatively, partly quan- 
titatively. We estimated the distance to the central en- 



14 Ziaeepour & et al. 



gine where the collisions between shells occurred to be 
~ 10 12 cm. This is consistent with the original suggestion 
of the model in which, internal shocks happen at a distance 
of r ~ 10 11 - 10 13 cm from the central engine. We also 
found that the external shocks influence the afterglow of 
this burst only at late times; most probably a few thou- 
sand seconds after the prompt emission. If our interpreta- 
tions are correct, it seems that the coherent magnetic field in 
the fireball can have relatively long lifetime and the energy 
distribution of electrons can be much more complex than a 
simple power-law. Based on the hardness ratio in the X-ray 
band we showed that there is a separation between popula- 
tion/sources of electrons responsible for the late X-ray and 
optical emission. In particular we conjectured the presence 
of a non-synchrotron, close to thermal component in the op- 
tical emission from the shocked heated material at late times 
in the afterglow of GRBs. The measurement of the late time 
optical spectrum can verify this hypothesis. A better under- 
standing of the relativistic plasma physics is also necessary 
and will provide realistic models for the evolution of fields 
and charged particles acceleration and energy distribution 
during the shock and their relation with the state of matter 
in the ejecta (fireball). 



APPENDIX A: 

In this appendix we briefly review the main aspects and 
results of the shock model used for comparison with data. 
It is a reformulation of the relativistic shock models with 
only one synchrotron emitting shocked layer but more de- 
tailed parametrization of the physical processes and their 
time variation (Ziaeepour, in prep.). 

We consider a radiative shock between two relativistic 
spherical shells. For an observer in the rest frame of the fast 
shell, the kinetic energy of the falling particles from the slow 
shell is instantaneously emitted as the shells join and merge 
together. A far observer sees that the fast shell absorbs the 
slow one, its density and/or size increases and its Lorentz 
factor decreases. During coalescence, we distinguish one re- 
gion of unshocked material in each shell, and one shocked 
zone in which charged particles are accelerated by the co- 
herent electric field produced be the shock. In presence of 
a coherent magnetic field, these charged particles irradiate 
their energy as synchrotron radiation. The magnetic field 
is usually conside r ed to be produced by the sh ock as well 
i|Yang et alj|l994 IWiersma fc Achterberd 120041 ) . It is also 
possible that there is an additional ambient magnetic field 
produced by the central engine. Here we neglect this possibil- 
ity. Thi s model is a simp lified version of the relativistic shock 
model (|Sari et alj|l996h that considers one shocked region 
in each side of the shock discontinuity - forward and reverse. 
The lack of a clear evidence for a reverse shock emission in 
any burst, specially during the prompt emission, means that 
the reverse shock in GRBs is weak and the assumption of 
just one shocked region is a good approximation. 

Energy-momentum conservation equations in the slow 



shell/ISM frame is: 
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where r is the distance from the central engine, n is the num- 
ber density of the fast shell measured in the slow shell frame, 
no is the number density of the slow shell in its rest frame, 
Ar is the thickness of the shocked synchrotron emitting re- 
gion, 7 is the Lorentz factor of the fast shell in the slow shell 
frame and j3 = \f"f 2 — I/7, m = m p + m e w m p , E 3y is the 
total emitted energy, and c is the speed of light. In the ab- 
sence of knowledge about the evolution of microphysics of 
the shocked region, we can not evaluate the evolution of Ar. 
However, the left hand side of equations I A 1 1 and I A2I depends 
on nAr - the column density of shocked region. 

The total power of synchrotron emission depends on the 
number of emitters. If only falling electrons are accelerated, 
the number of emitting sources are proportional to no- In 
this case we have to consider a model for the evolution of 
the shocked region. A minimal estimation is Ar = /3Aro/j 
where Aro is a constant thickness scale. The physical rea- 
son for this choice is that in the two extreme cases when 
the relative Lorentz factor is very large or very small, one 
expects a thin shocked region. In the first case the reason is 
that due to the large density difference between shells only a 
thin layer of the fast shell is affected (shocked) be the falling 
particles from slow shell. In the second case the shock is soft, 
the two shells merge smoothly and one expects that the size 
of the turbulent shocked particles be very restricted. An- 
other plausible model is the acceleration of all the electrons 
in the shocked region following the formation of a coherent 
electric field in the shock. In this case, the column density 
of emitters is proportional to nAr, and therefore an ad- 
hoc expression for the evolution of Ar is not needed. Both 
these cases are based on the assumptions that only a detail 
knowledge about the micro-physics of the matter state and 
shock can confirm or role out. Nonetheless, if the rest frame 
densities of the shells are close to each other, these approx- 
imations should give similar results. For simplicity, here we 
will consider the first case which makes the right hand side 
of I A II and I A2I independent of n. 

We expect that Ar/r <C 1 although in some circum- 
stances the thickness of the shells can be comparable with 
r. If r — > 00, the shells become planar. If the shells are not 
spherical but collimated, the formulation of the shock and 
synchrotron radiation is more complex. However, it can be 
shown that at first approximation, the total flux can be con- 
sidered as proportional to the opening angle of the jet. Hence 
for simplicity we do not discuss the effect of the collimation 
further here. 

Variation of the synchrotron energy can be related to 
the emission power dE sy /dr — dE ay //3dt = P sy /f3, where 
P 3y is the total power of the synchrotron emission. Assum- 
ing a power-law distribution for the Lorentz factor of the 
accelerated electrons, dE ay /dr can be written as a function 
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of the relative Lorentz factor 7, e e , es, and Jj], respectively 
the fraction of kinetic energy of the shell transferred to accel- 
erated electrons, the fraction of kinetic energy transferred to 
a coherent magnetic field, and the number distribution index 
of accelerated electrons l|Sari et al.lll998l ). We assume that 
q is a constant but e e , es have power-law dependence on r 
with indices a e and olb respectively. Then we can solve equa- 
tions IA1I and IA2I to find the evolution of 7 and the column 
density of synchrotron emitting layer nAr. These equations 
are coupled and therefore some approximations are neces- 
sary to decouple and solve them. In addition, their solutions 
depend on the relative strength of the various parameters. 

At lowest order in e, and when r/ro — 1 <C 1, the evo- 
lution of 7 can be expressed in the following form: 



7 2 W 
7 2 (r ) 



1 + Fe + ... ro 



The coefficient F depends on various physical quantities 
such as initial relative Lorentz factor of the colliding shells, 
their initial densities, e e , eb, the number distribution index 
of accelerated electrons q, and the variation of these quan- 
tities with time. They are degenerate and the extraction of 
all of them from available data is not possible. 

The column density of synchrotron emitting layer at the 
lowest order and with approximations mentioned above is: 
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where the term r%jr 2 is simply due to the adiabatic expan- 
sion of the shell. The constant F\ depends on various param- 
eters mentioned about. The exponent a 1 consists of 2 parts: 
a constant and a part dependent on r\ = 2a e + as. This 
shows the importance of the evolution of coherent electric 
and magnetic fields in the GRBs. 

Evidently, the only observable outcome of the shell col- 
lision for us is the synchrotron emission. Its in tensity for one 
charged particle is well known l|jacksonlll999h . and in order 
to obtain the total intensity we should integrate it over the 
emitting volume and the spectrum of the synchrotron emit- 
ting charged particles (electrons): 



dl r- e 2 f°° 

— = 2V37T — V(r)7 / d7 e n e (7 e ) 

av c J r 



Ks{x)dx (A5) 



where V(r) is the volume of the synchrotron emitting layer, 
v c is synchrotron characteristic frequency, n e (7 e ) is the elec- 
tron density with Lorentz factor 7 e , Y is the minimum 
Lorentz factor of the accelerated electrons, and finally K 5 

is the Bessel function. Although 4£ is the main observable 
quantity, it depends on too many unknown parameters and 
in practice it can not be efficiently used. Two other quanti- 
ties derived from intensity are more useful. One is the hard- 
ness ratio and the other is the lag of the peak emission in 



7 The more popular parameter p the energy distribution index of 
the accelerated electrons is related to q: p = q — 1. 



different energy bands: 



HRi 



Ae, ; 



peak 



peak 



d_ 

7k 







(A6) 



(A7) 



where i and j are different energy bands with logarithmic 
mean frequency v and minimum and maximum frequency 



and v" 



respectively. The variable e is defined in 



IA3I The hardness ratio can be determined directly from the 
data for all energy bands. However, the lags are only mea- 
surable for fast varying features such gamma-ray peaks and 
X-ray flares. Nonetheless, when they are measurable they 
give valuable information about the distance to the central 
engine at which, the corresponding events have occurred. 
Note that e's in IA7I are dimensionless relative quantities. 
However, the measured lags are not. Therefore, if the lags 
in multiple bands are available, we can estimate the value 
of the corresponding e's thereby we can measure the initial 
distance ro- 

It can be shown that HRij and Ae^ are depend on r (or 
equivalently t or e) through the function C(r, v) = v/v' c Y 2 — 
v/vm, where v' c = eB/m e c depends only on the magnetic 
field. Therefore, even without any knowledge about other 
quantities, the evolution of the hardness ratio directly shows 
the evolution of v m . The analytical expression for HRij and 
specially for Aeij are quite involved. Moreover, the integral 
over the Bessel function in IA5I must be calculated numeri- 
cally. Our experience shows that the use of asymptotic val- 
ues of the Bessel function which permit to obtain analytical 
expressions for these quantities leads to large errors and un- 
physical behaviour of the results. The results of numerical 
calculation for a number of combination of parameters will 
be reported in (Ziaeepour, in prep.) and we have used them 
for analyzing GRB 060607A data. 



APPENDIX B: 

Practically in all gamma-ray bursts the spectrum becomes 
softer after the end of the main prompt spikes. Therefore 
the spectral index of the gamma-ray spectrum should be 
considered as a slowly time-varying quantity. The total X- 
ray flux / 7 _x(t) expected is approximately: 



fi-x{t) 



(Ex 

f (t) ( e \ a T (t) 

J1\") I -^Tnax \ 



dE* 



a~,(t) + 1 \ E, 



(Bl) 



where is the logarithmic mean of gamma-ray energy 
band, E max and E m i n are maximum and minimum of the 
X-ray energy band, and a 7 (t) is the slowly varying spectral 
index. 

In the method explained in Sec.|3]for the extrapolation 
of BAT spectrum to X-ray band, a simple average of the 
spectral index, usually from the total BAT spectrum and 
early X-ray spectrum is used as an approximation to take 
into account the time evolution of the spectrum. Using the 
same formulation as in equation IB1I the extrapolated BAT 
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light curve to the XRT band fBAT-XRr(t) can be estimated 
as: 

t u\ ~ Nf B AT(t) ( E max \ a 

j B AT— XRT (t) « m + 1 \~- ) (B2) 

where a is the average of BAT and XRT spectrum indices 
and N is a normalisation factor between the two instru- 
ments. In general, the cruder approximation of IB2I will not 
agree with the more precise expression IB 1 1 unless a 1 {t) have 
roughly a linear and slow evolution with time. If these con- 
ditions are not fulfilled and/or the prompt and the tail emis- 
sion in X-ray have different origins, the extrapolated light 
curve lB2l should deviate from the observed one, or a break 
appears when one tries to join simultaneous observations by 
BAT and XRT. This break can not be removed by adjust- 
ing the nor malisation factor N in lB2l Among 40 Swift GRBs 
studied in l)0'Brien et al.ll2006h only a few of them seems to 
need such a break to join the extrapolated BAT and the ob- 
served XRT light curves: GRB 050315, GRB 050713B, and 
GRB 050915B. In conclusion. IB2I is a good approximation 
of IB II for most bursts and when it is not, the probability of 
having a smooth and overlapping light curve with the XRT 
observations is very small. 
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